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Abstract 

Rates for nucleation of two-flavor quark matter in a neutron star core, originally 
composed of nuclear matter in P equilibrium, are calculated at zero temperature by 
a quantum tunneling analysis incorporating the electrostatic energy. We find that a 
nucleated droplet would develop into bulk matter due to electron screening effects. 
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The possibility that quark matter could exist in neutron stars and the astrophysical 
consequences have been considered for the past two decades. (See, e.g., Refs. l)-3).) Usually, 
the main question was at what pressure the free energy per baryon for electrically neutral 
quark matter becomes less than that for nuclear matter. Recently, Glendenning^^ discovered 
that, by relaxing the constraint of local charge neutrality, a phase where quark and nuclear 
matter in f3 equilibrium coexist in a uniform sea of electrons could appear for a finite range 
of pressures. This is because the presence of strange and down quarks plays a role in 
reducing the electron Fermi energy and in increasing the proton fraction of nuclear matter. 
His work^) and the subsequent work by Heiselberg et al.^^ claimed that such a mixed phase 
could exhibit spatial structure such as quark matter droplets embedded in nuclear matter 
due to surface and Coulomb effects. As the star whose core consists of nuclear matter in 
f3 equilibrium spins down or accretes matter from its companion star, the central density 
could become sufficiently large for the mixed phase to be stable. Whether the mixed phase 
actually nucleates, however, depends on the occurrence of the dynamical processes leading 
to its nucleation. The first of these processes should be the conversion of nuclear matter to 
two-flavor quark matter, since the conversion to three-flavor quark matter, which requires 
many simultaneous weak interactions, is unlikely to occur. In this paper, therefore, we 
consider at what pressure a droplet of two-flavor quark matter forms in nuclear matter and 
whether it develops into bulk matter or remains finite. 

We begin by describing the bulk properties of the various components for densities well 
above the nuclear saturation density. For nucleons, we write the simple formula for the 
energy density adopted by Heiselberg et al.:^^ 



where m„ is the neutron mass, n is the nucleon number density, Uq — 0.16 fm~^ is the 
nuclear saturation density, x is the proton fraction, Kq — 250 MeV is the coefficient of the 
compressional term, and Sq — 30 MeV as well as 7 = 1 determine the symmetry term. For 
quarks we adopt the energy density based on the bag model. 



where q — u,d, and s denote up, down, and strange quarks, Uq is the number density of q 
quarks, and jj,q = {1 — 2Q;s/7r)~^/^(7r^nq)^/^ is the chemical potential of q quarks. In Eq. (2), 
all quark masses have been taken to be zero, and the QCD fine structure constant and the 
bag constant have been set as — 0.4 and B — 120 MeV fm~^. The energy density of the 
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electrons, which are relativistically degenerate, is ee(ne) = fif/An'^ with the electron number 
density Ue and the electron chemical potential fig = (STT^ng)^^'^. We neglect the existence of 
hyperons and muons. 

At a given pressure P, we have calculated 
the baryon chemical potentials for the three bulk 
phases being electrically neutral. The first phase 
consists of nuclear matter in f3 equilibrium: /i„ = 
l^p + /^e with the neutron (proton) chemical poten- 
tial Hn ifJ'p)- The second phase is composed of u and 
d quark matter arising via deconfinement from the 
above nuclear matter of proton fraction 
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Lastly, for com- 



parison, we consider a phase containing u, d, and 
s quark matter in (3 equilibrium: fi^ = fJ'u + fJ'e and 
/^rf = /^s- Figure 1(a) shows the chemical potentials 
for these three phases obtained as jj, = {e + P) / rii, 
with the corresponding total energy densities e and 
baryon densities ni,. The crossings denote the tran- 
sition points from one bulk phase to another, at 
which the baryon density jumps as can be seen from 
Fig. 1(b). 

Let us now inquire how long it takes a u and d 
quark matter droplet to form via deconfinement in 
/3-equilibrium nuclear matter for pressures in the 
vicinity of the transition point shown in Fig. 1(a). 
We assume such a droplet to be a sphere macro- 
scopically characterized by its radius R and de- 
scribe the tunneling behavior of a virtual droplet 
in the semiclassical approximation. A Lagrangian 
for the fluctuating droplet may thus be written as^^ 

1 




Fig. 1. (a) Chemical potentials for the 
electrically neutral bulk phases as 
a function of pressure. The solid 
line stands for the phase of nuclear 
matter in P equilibrium; the dashed 
line, the phase of u and d quark 
matter arising therefrom via decon- 
finement; the dotted line, the phase 
of u, d, and s quark matter in (3 
equilibrium, (b) Energy densities 
for the corresponding phases as a 
function of baryon density. The 
double-tangent constructions (dot- 
dash lines) denote the coexistence 
of the two bulk phases involved. 
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where M{R) is the effective droplet mass, and U {R) is the potential for droplet formation. 
In the present analysis, no dissipation is considered. 

We proceed to evaluate the potential U{R). Due to the high sound velocity of the system 
(~ c) it may be assumed that the number density of each component adjusts adiabatically to 
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fluctuations of R. Consequently, the system retains pressure equilibrium between quark and 
nuclear matter, and each component, if it makes no contribution to Coulomb screening on 
the resulting excess droplet charge, is homogeneous, or else is distributed accordingly. Here 
we assume the quark distributions to be uniform; their shifts due to quark screening^) should 
reduce the Coulomb energy only slightly. The electron gas is most effective at screening; 
its Thomas-Fermi (TF) screening length is ~ 6 fm since ~ 300 MeV. For the droplet 
sizes of interest i? < 10 fm, protons, whose TF screening length is estimated from Eq. (1) 
as > 10 fm, may be taken to be uniform. For the initial metastable phase of nuclear matter 
in j3 equilibrium under pressure P, we set the chemical potential as yUjuit = (ejnit + -P)/'^init 
with the energy density ei^it = tN{ni^\uX;,,iO +^e{ne,\nit), where ni„it, Xjnit, and rieinit give the 
initial number densities of nucleons and electrons. We can thus express the energy density 
for the inhomogeneous phase containing a single droplet as the sum of bulk, surface, and 
Coulomb terms: 

en = 0{R - r)tQ{nu, na, 0) + 9{r - -R)eAr(ninit, Xinit) + ee[ne(r)] + es(r) + — - — , (4) 

OTT 

where r is the distance from the center of the droplet, es{r) is the increase in energy density 
due to the quark and nucleon distributions in the surface layer whose thickness we assume 
to be much smaller than R, and E{r) is the electric field. We then integrate the difference 
in the thermodynamic potential per unit volume at chemical potential //init between the 
initial and inhomogeneous phases over the system volume V which is related to P as P = 
-d{Vei^it)/dV. We thus obtain U(R) as 

AttR^ 

U{R) = -^nb,Q{fiQ - /iinit) + AEe{R) + ATTaR"" + Ec{R) , (5) 

where — (n„ + nd)/S is the baryon density inside the droplet, fj,Q = [eQ(n„,nd) + 
ee(ne,init) + IS the chemical potential for the uniform quark matter determined by n„, 

nd, 'T'e.init, and P, AEe{R) = Jy dV{ee[ne{r)] — ee('T'e,init)} is the excess of electron energy over 
the initial one, a is the surface tension, and Ec{R) is the electrostatic energy. In Eq. (5) 
the curvature energy (oc R) is omitted. Estimates of the quantities n„, n^, ne{r), AEg{R), 
(7, and Ec{R) are as follows. 

The distribution of electrons has been determined in the linear TF approximation. Elec- 
trons act to screen the excess droplet charge Ze = (47rP^/3)pQ with pq = e[(2n„ — nd)/3 — 
^e,imt], and thereby deviate from uniformity in proportion to the electrostatic potential 0(r) 
as Sne{r) = ne(r-) — ne,mit = {due/ dfie)imte(p{r) . The potential 0(r) satisfies the Poisson 
equation V^0(r) — ^^(^{r) = —A'kpq9{R — r), where n — J ATie^{dne/ diJie)imi is the inverse of 
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the TF screening length. Its solution reads as 

(47rpQ/«;^)[l — exp(— + kR) sinh(«;r)/Kr] , 



r<R 



(f){r) 



(6) 



{Attpq/ K^)[KRcosh{KR) — smh{KR)] exp{—Kr)/Kr , r > R . 



Figure 2 exhibits a typical example of the electron screening at i? ~ 3 fm. When the initial 
chemical potential of electrons //e.init ^ e0, as in Fig. 2, the linear TF approximation can be 
safely used. For i? < 10 fm, it turns out that e0///e,init ^ 0.1. Within the confines of this 
approximation, AEg{R) — Z^e,imt- 

We then calculate the electrostatic energy Ec{R) using Eq. (6). The result is 

Ec{R) = (27rV|//«^){-3+ («;i?)^ + exp(-2«;i?)[3 + 6Ki? + 5(«;i?)2 + 2(«;i?)2]} . (7) 



Note that, in the limit of electron uniformity {kR - 
The surface tension cr, being poorly known, has 
been taken from the Fermi-gas model for a quark 
matter droplet in vacuum.^) Due to the reduction in 
quark density of states arising from the droplet sur- 
face, the surface tension a — (3/47r^)(TOu/x^-|-md//^) 
occurs for ^ rUq {q — u,d). Since ~ 500 
MeV, Hd ~ 600 MeV, and niu ^ ma ^ 10 MeV, 
the value of a is roughly estimated as cr ~ 10 MeV 
fm~^. The Fermi- gas description of finite quark 
matter as adopted here, however, is problematic 
for a small baryon number inside the bag ( < 5).^^ 
Nevertheless, we expect such a description to be 
valid, since the baryon number of a virtual droplet 
moving under the potential barrier is > 10 (typi- 
cally ~ 100). 

Finally, we determine the number densities 
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Fig. 2. Deviation Sn^ of the electron 
number density from the uniform 
one ne,init due to the TF screening 
on the excess droplet charge of den- 
sity pQ and of radius R aX P = 455 
MeV fm~^. See the text for the def- 
inition of r, K, //e,mit) and (f). 



and rid- As a result of deconfinement, these are re- 
lated via riu/rid — (l+a^init)/(2 — a^^init)- Another relation comes from the pressure equilibrium 
between quark and nuclear matter, 

AE,{R) 



N 



AttR^ 



(8) 



where Pq = [eQ{nu,nd) — 4i?]/3 is the quark pressure, and Pjv — P — ee(ne,init)/3 is the 
nucleon pressure. As long as cr oc n^/^, the surface tension contributes nothing to the 
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pressure equilibrium (8). In Eq. (8) the Coulomb pressure arising from Ec{R) has been 
ignored since it is negligibly small. n„ and nd are thus independent of R. 

We next estimate the effective mass M{R) from the kinetic energies of nucleons and of 
electrons. We obtain the nucleon kinetic energy by solving the nucleon continuity equation 
with the boundary condition at the droplet surface^) as 27reiv(^mit5 3;init)(l — '^6.Q/'"'init)^-R^-R^- 
The electron kinetic energy, which is calculated from the electron continuity equation as 
< 27ree(ne,init)(pQ/'"'e,mite)^i?'^-R^, proves negligible since e^v ~ 1000 MeV fm^^ 3> eg ~ 40 
MeV fm~^. We thus obtain 

MiR) = 47rejv(ninit, a^imt) (l - ^) ' R' ■ (9) 

V ?^init/ 

We turn to the evaluations of the time r required to form a single droplet. By calculating 
not only the energy Eq for the zeroth bound state around R = but also the corresponding 
oscillation frequency and probability of barrier penetration po from the Lagrangian (3) 
in the WKB approximation the values of r have been obtained as r = {i^oPo)~^- With 
increasing pressure from the transition point Pq = 445 MeV fm~^, the potential barrier 
is lowered for fixed a as the sohd hnes in Fig. 3 exhibit. Since this barrier controls the 
underbarrier action Sq between the classical turning points and hence the probability po ~ 
exp{—2So/h), the nucleation time r calculated for a — 5, 10, 15, and 20 MeV fm~^ shows an 
exponential P dependence in the overpressure regime, as illustrated in Fig. 4. The oscillation 
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Fig. 3. Potentials for the formation of a 
droplet with the surface tension cr = 10 

MeV fm-2 at P = 445, 450, 455, and 460 
McV fm^'^. The sohd lines are the results 
for the electron-screened droplet and the 
dashed lines are those for the non-screened 
droplet. 




440 450 460 470 480 
P (MeV fm-3) 



Fig. 4. Time for the formation of a single 
droplet with the surface tension a = 
5, 10, 15, and 20 MeV fm~^ as a function of 
pressure. Pq denotes the transition point. 
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time scale is naturally ~ 10^^"^ sec. The energy Eq is mostly comparable to the barrier 
height ~ 100 MeV. Nevertheless, the WKB approximation is expected to be useful because 
the WKB analysis of tunneling for the Lagrangian having no Coulomb energy but otherwise 
the same R dependence as Eq. (3) agrees well with the fully quantum analysis. It is 
important to note that the formation time of the first droplet anywhere in a neutron star 
is t/N, where is the number of virtual centers of droplet formation in the star. N is 
uncertain, but is at most the total baryon number in the star (typically ~ 10^''). 

Finally, we give some astrophysical implications. The nuclear equation of state of inter- 
mediate stiffness used here allows a neutron star to have the central density higher than the 
transition point ~ Sno.^^-* Once the central pressure reaches the pressure giving a realistic 
value of t/N during the spin-down or accretion, a droplet of u and d quarks may appear in 
the nuclear matter. We can see from Fig. 3 that the electron screening, without which Eq 
(oc R^) dominates U{R) for large R, prevents the droplet from remaining finite for i? < 10 
fm. For R ^ Kr^, Eq behaves as Eq oc R^ and Eq < 'inaR^^ hence U{R) < 0. The droplet 
may thus expand into bulk matter with the speed of sound, induced by the density jump 
due to deconfinement. The resultant bulk quark matter occupying the region P ^ Po may 
change into three-fiavor quark matter via weak interactions in a time scale of 10~^ sec.^^^ As 
long as three-flavor quark matter is more favorable than nuclear matter, the region of u, d, 
and s quarks may spread^^^ with a release of chemical energy which amounts to ~ 100 MeV 
per baryon as estimated from Fig. 1(a) and turns immediately into thermal energy. This 
implies that the mixed phase consisting of quark and nuclear matter^''^' is unlikely to occur. 
The resulting star could be stable, so the thermal energy deposited in the star ~ 10^^ ergs 
might be a possible origin of a 7-ray burst. 

In summary, we have found that a droplet of u and d, quarks, if it appeared in a neutron 
star, would develop into bulk matter due to the electron screening effects. The effect of 
muons, which was neglected above, makes nuclear matter more stable (Pq ^ 465 MeV 
fm^''^) and shortens the TF screening length by ~ 1 fm. We have confirmed that such 
a change makes no important difference. In order to make better estimates, we should 
examine dissipation effects on the nucleation, which may depend on the superfluidity of 
nuclear matter (e.g., Rcf. 16)) and control the crossover temperature from classical^'''' to 
quantum nucleation, the roles of hyperons in the composition of matter before and after 
deconfinement, and the poorly- known interfacial properties, the effects of curvature from 
the quark side being shown to significantly destabifize the droplet. In addition to recent 
low-temperature experiments,^^) astrophysical dense matter such as neutron star matter of 
temperature < 0.1 MeV may give us examples of quantum nucleation at first-order phase 
transitions. 
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